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Abstract 

We have proposed a qualitative model for the structure of binary systems similar to Pd-Er 
alloys, which explains their nonmonotonic relaxation after the hydrogen saturation. It is based 
on the assumption that such a solid solution involves two kind heterogeneities. The former 
are caused by spinodal decomposition of the initially homogeneous state of the solid solution 
into the phases enriched and depleted of Er atoms. The latter are crystalline defects that trap 
an additional amount of Er atoms, which leads also to their local accumulation, changing the 
defect properties. Hydrogen atoms penetrating into the solid disturb the equilibrium of both 
the phase separation and the defect saturation with Er atoms, causing redistribution of Er 
atoms. The diffusion fluxes give rise to the motion of the interface between the two phases 
that is responsible for time variations, for example, in the relative volume of the enriched phase 
observed experimentally. We have found the conditions when the interface motion can change 
the direction during the system relaxation to a new equilibrium state. The latter effect is, from 
our point of view, the essence of the hydrogen induced nonmonotonic relaxation observed in 
such systems. The numerical simulation confirms the basic assumptions. 



1 Introduction 

Metal-hydrogen systems are singled out with the fact that they can remain in nonequilibrium states 
for a long time because of the ability of hydrogen atoms to migrate relatively fast through metallic 
matrices. Therefore, a metallic matrix initially saturated with hydrogen inevitably loses its main 
amount during a certain time period which may be sufficiently long fl], ^ §]• Hydrogen atoms 
penetrating into the metallic matrix cause the reformation of the defects structure and as a result, 
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the system structure becomes substantially heterogeneous, moreover, it can be of multiscale geometry 
^. When the heterogeneity scales expand over several levels these systems may be treated as 
fractals (see, e.g., 0, 0; i|) and their evolution goes in a sophisticated manner. In particular, under 
certain conditions the dynamics of such a systems is described by strange attractor, which reflects 
in nonmonotonic or even chaotic time variations in its structure and occurring phase transitions |T0|] . 

The characteristics of structure transformations during the hydrogen saturation and the following 
relaxation have been studied for metallic alloys based on Pd fl^, |lT], [l^, |13], [l^. For example, for 



the Pd-Er alloys undergone deformation the hydrogen saturation gives rise to the change in the sign 



of the elastic tension at the initial stage of relaxation This effect is due to variations in the 

image forces that are induced by the transformation of the defect-metal (D-M) complexes into the 
hydrogen-defect-metal (H-D-M) ones because of the high binding energy between hydrogen atoms 
and defects in Pd matrix |T^. Having changed the sign, the tensions grow during two days, then, 
drop and decrease by 25% in eight days. In parallel with it oscillatory variations in the relative 
volume of the phases enriched and depleted of Er atoms occur practically right after the hydrogen 
saturation. The difference in the concentration of Er atoms between the two phases evolves in the 



same manner. These oscillations are chaotic rather than regular |Ty, |T^ . 

The given phenomena were observed in such Pd-Er alloys where the atomic concentration of Er 



atoms was about 10% [T^. So their structure has to be sufficiently simple, namely, should comprise 
only two phases formed of the Pd matrix depleted of or enriched with Er atoms. So the relaxation 
after the hydrogen saturation seems to go only through time variations in the relative volume of 
these phases without complex structural transformations in the Pd-Er system that should appear 



for high values of the Er concentrations |[T6 



The stochastic behavior of the relaxation process in the Pd-Er-H system has been demonstrated 
in a frameworks similar to the Lorenz model proposed for analysis of turbulent phenomena in atmo- 
sphere IjT^ and, then, applied to the description of plastic deformation [|T^. It should be pointed 
out that time variations in the structure of the given alloys after the hydrogen saturation enable us 
to consider diffusion flux in them also turbulent. It can be justified at least for some stages of the 
relaxation if we take into account not only the stochastic evolution of the alloy structure but also 
the difference in the rate of nonmonotonic structure transformations for different coherent diffraction 
regions, as it was noted for the first time in |TB]. The analysis presented in |TD] has been carried out 



in the frameworks of macroscopic approach, so, in order to complete it, we should develop a micro- 
scopic mechanism by which these phenomena can be implemented. The latter problem is actually 
the purpose of the present paper. 

Papers |jlO|, |13| demonstrate that the macroscopic model to be developed for the hydrogen induced 
relaxation in alloys similar to Pd-Er-H has to take into account the following. 

1. In such systems the hydrogen saturation initiates several interacting processes, in particular, the 
transformation of the D-M complexes into D-M-H ones and the diffusion of Er atoms between 
the enriched and depleted phases of the Pd-matrix. These two processes cause nonmonotonic 
time variations not only in the relative volume of these phases but also in the Er concentration 
inside them. 

2. The correlated variations in the Er concentration and the phase relative volume suggest that 
the two processes stem from the same phenomenon. Besides, their nonmonotonicity points out 
that diffusion of Er atoms does not obey conventional Pick's low, at least, everywhere in the 
alloy. 
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3. The D-Er-H complexes induced by the hydrogen saturation have to be enriched not only with 
hydrogen atoms but also with Er atoms because of high affinity of H and Er atoms. So we 
meet four type regions (fig. 2): two regular phases of the Pd- matrix and regions with two 
type defects D-Er-H (inside each of the phases). Besides, substantial difference in the specific 
volumes of the regular phases and complexes can lead to the appearance of transitional regions 
between the matrix and complexes. 

4. The complexes are the source of additional nonuniformity in the distribution of Er atoms, 
causing an additional diffusion flux in the system. 

5. The variations in the phase relative volume mean displacement of the phase boundaries among 
with diffusion of Er atoms. So dealing with the nonmonotonic relaxation on systems similar to 
Pd-Er alloys we first of all should explain why a nonmonotonic motion of the phase boundaries 
occur. 



It should be pointed out that the system Pd-Er-H is not unique keeping in mined the aforesaid 
properties. Even more pronounced chaotic behavior of the structural evolution has been observed 
in the Pa-Ta-H system [0, which suggests the multiscale organization of its structure and the 



substantially nonequilibrium state of the defects and the adjacent regions of the Pd-matrix. 

In spite of the available number of theoretical models cited above, physics of the hydrogen induced 
relaxation in such metallic systems is far from being understood well. So at the ffist step of developing 
a microscopic model able to explain the aforesaid phenomena, it is reasonable to confine ourselves 
to the analysis of the interaction between the D-M-H complexes and the phase boundaries through 
diffusion of Er atoms. In this study we should also take into account that the distribution of Er atoms 
around the D-M-H defect can be nonmonotonic because under certain conditions the deformation 
field in their neighborhood exhibits a nonmonotonic behavior. Nevertheless in the present paper we 
consider only defects with a sufficiently simple structure and propose a model which qualitatively 
explains the nonmonotonic motion of the phase boundaries due to the hydrogen saturation. We do 
not claim on a sophisticated theory of this relaxation process but only show the feasibility of its 
qualitative explanation in the frameworks of the proposed model and single out the main clues to its 
understanding. 



2 Model 



The model to be developed uses the basis of the molecular-kinetic theory of solid solutions ^ ^2 
^ and the state diagram of the Pd-Er alloys (in the cT-plane, where c is the atomic concentration 
of Er atoms, T is temperature) |jl6[. Namely, we assume that difference in the atomic radii and 
electro-chemical characteristics of Er and Pd results in instability of the homogeneous state of the 
given solid solution Pd-Er. The instability gives rise to domains Qet with an increased value c+(T) 
of the Er concentration whose composition approximately corresponds to the intermetallic ErPd7, 
the ffist equilibrium phase in the state diagram after the primary solid solution of Er atoms in the 
Pd-matrix |[T^. These domains are surrounded by the maternal phase Qpd (primary solid solution) 
with a lower value c_(T) of the Er concentration equal to the limit of Er solubility in the Pd-matrix. 
In equilibrium for the regular solid solution the relationship between the values c+(T) and c_(T) is 



specified by the "arms rule" pOf. However crystalline defects such as vacancies, interstitial atoms. 
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their complexes, dislocations, grain boundaries, and so on can trap Er atoms. Furthermore we 
cannot exclude that other phases more enriched with Er atoms, ErPds, ErioPd2i, etc., appear in the 
neighborhood of these defects. A multiscale hierarchical structure can form itself in this system in 
general. 

Penetrated into the alloy hydrogen atoms disturb the equilibrium in two ways. First, the hydrogen 
saturation converting the solid solution Pd-Er into the ternary system Pd-Er-H changes the phase 
equilibrium in the metal subsystem. The higher affinity between H and Er in comparison with the Pd 
and H should give rise to the complexes {Er-H} that contain not only binary bonds Er-H but also 
many bond structures, for example, Er-H-Er. The first type clusters, Er-H, seem to be typical for 
regions with low values of the Er concentration, whereas in the Er enriched phase Qsr as well as near 
the solubility limit of Er atoms in the Pd-matrix the amount of the second type clusters should be 
essential. The presence of the Er-H-Er bonds effectively makes the interaction of Er atoms with each 
other stronger, which has to cause the state diagram of the triple solution Pd-Er-H take a form shown 
at fig. 1^. Drawing this diagram we were keeping in mind one presented in [|I^] for the binary solution 
Pd-Er and treating hydrogen atoms as a source of the additional Er-Er attraction. Let us explain 
the proposed form of this diagram in more details. Under the adopted assumptions the instability 
of the system homogeneous state is to occur at lower concentrations of Er atom than it would be 
in the absence of hydrogen atoms. Correspondingly, the solubility limit c_(ch,T) of Er atoms in 
the Pd-matrix has to depend on the hydrogen concentration decreasing with the value of ch grows. 
In its turn, the composition of the intermetallic ErPd2^(H) or, what is the same, the value c+(ch,T) 
seems also to depend on the local values of the hydrogen concentration ch going into the region of 
higher value with ch grows. We assume the intermetallic ErPda;(H) to originate from ErPd7 and to 
form a certain phase in the nonequilibrium system under consideration. We note that the concept 
of true intermetallic concerns a rigorously ordered phase with specific stoichiometric composition. 
So the concept of intermetallic for nonequilibrium systems requires an individual consideration. For 
given system the experimental data demonstrate the presence of a certain phase containing hydrogen 
atoms where the concentration of Er atoms closely matches the composition of the intermetallic 
ErPdy. However, whether this phase can be treated as a true intermetallic ErPdy or it comprises 
its small clusters, equilibrium or nonequilibrium, containing hydrogen atoms is of low importance 
in the present study. We need only the presense of a certain quasiequilibrium phase experimentally 
observed, where the atomic concentration of Er atoms is about 10%. We suggest this phase to be 
relative to the true intermetallic ErPdy so below we also will refer to it as to "intermetallic" enclosing 
with quotation marks. 

Second, due to higher affinity between hydrogen and the defects in comparison with the regular 
structure of the Pd-matrix hydrogen atoms have to be attracted by defects containing also an in- 
creased amount of Er atoms. By the same reasons as discussed in the previous paragraph hydrogen 
atoms penetrating into the defects located in the phase Qpd have to density the Er distribution in 
them, and, so, these defects gain the ability to absorb new Er atoms. In another words, hydrogen 
atoms penetrating into phase Qpd increase the "capacity" of defects being activated by them as traps 
or sinks of Er atoms. Hydrogen atoms penetrating into the "intermetallic" Qet also affect the state 
of its defects make variations in the "intermetallic" composition possible. In fact, without hydrogen 
the composition of ErPdy is fixed (in the state diagram the intermetallic ErPdy corresponds to a 
line rather than two-dimensional region [T^). So, defects located in this intermetallic are to contain 
other phases enriched with Er atoms to higher degrees, e.g., ErPds. Now, these local phases gain the 
possibility to dissolve in the "intermetallic" Q-^r injecting new Er atoms into in it. In other words. 



4 



the hydrogen saturation into the phase Qsr activates its defects as local sources of Er atoms. 

The interface T separating the phases Qpd and Qei is of atomic size. Therefore during the 
relaxation course this interface locally must be in quasiequilibrium, i.e. the concentration of Er 
atoms in the regions of the phases Qsr and Qpd adjacent to the interface T must be equal to 
c+(ch, T) and c_(ch, T), respectively. This quasiequilibrium is the result of practically instantaneous 
redistribution of Er atoms on atomic scales in the vicinity of the interface T. The following diffusion 
of Er atoms into the bulk of the phases Qev and Qpd tends to make their individual composition 
uniform, causing the interface T to move. In fact, generally the diffusion flux at the interface T on 
the side of the phase Qpd and that of the phase Qsr is different. Thus, due to atom conservation the 
interface T is to move to compensate this difference giving rise to an effective source of Er atoms 
proportional to f [c+(ch,T) — c_(ch,T)], where v is the normal velocity of the interface T. 

Nonmonotonic motion of the interface T being the essence of the observed relaxation process is 
due to change in the profile of Er distribution as time goes on. Fig. ^ sketches the Er distribution at 
the initial stage of relaxation after the hydrogen saturation. Light arrows point the diffusion direction 
of Er atoms in phases Qpd and Q^r- Whence it follows that for each of the phases individually the 
change in the Er equilibrium concentration near the interface T and the presence of the activated 
defects cause the diffusion flux in the opposite directions. At the initial stage the distant defects 
cannot affect the interface motion, so, it is governed by the change in the equilibrium state only. 
As time goes on the diffusion flux induced by the activated defects can become dominant and the 
interface T will move in the opposite direction. 

We suppose that this phenomenon is a backbone of the nonmonotonic relaxation induce by the 
hydrogen saturation into systems like the Er-Pd alloy. 

Before passing to a specific mathematical model we would like to note the following. As seen 
from the present analysis the individual role of the two phases is approximately the same. The latter 
enables us to confine our consideration to only one of them, for example, the phase Qpd- Taking into 
account both of them only numerical coefficients in the expressions to be obtained can be changed, 
makes no sense at the current stage of the theory development. 



3 Master equations 

Taking into account aforesaid we write the diffusion equation for the atomic concentration c of Er 
atoms in the phase Qpd regarded as the half-space z > (we have attached the coordinate system 
to the interface T moving at the velocity v in the laboratory frame): 

dc dc d^c 1 

where D is the diffusion coefficient of Er atoms treated as a constant and the second term in the 
right-hand side of equation (|^) describes the Er atom trapping by a defect approximated by the 
5-function located at a point ztr- This trap-defect is characterized by the following parameters, / is 
its real physical size, q is the capacity, i.e. the number of free seats for Er atoms at the current time, 
and Ttr is the characteristic time during which, on the average, the defect traps an Er atom located 
in its neighborhood (not to be confused with the lifetime of Er atoms inside the defects regarded 
infinitely long in the given model). In the frame {z,t} attached to the interface T the trap-defect 
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moves towards the interface at the velocity v, i.e. 

(2) 



dZtr 



dt 

Here z^j. is the initial position of the defect (at t = 0) and is one of the model parameters. Since the 
newly trapped atoms occupy vacant places inside the defect its capacity decreases in time according 
to the equation: 

dq _ qlc{ztr,t) 
dt aTti 

and go; the initial capacity of the trap-defect activated by the hydrogen saturation, is also a model 
parameter (a is a period of the lattice). 

At the interface T {z = 0) the Er concentration is set equal to: 

Ci = Co - Apd , (4) 

where cq is the initial concentration of Er atoms in the phase Qpd before the hydrogen saturation and 
Apd is a change in this concentration caused by the hydrogen effect on thermodynamical equilibrium 
of the phase boundary. Ignoring diffusion of Er atoms inside the phase Qei we can represent the law 
of Er atom conservation at the interface T in the form 



^ dc 
Av = D — 

oz 



(5) 



where A = c+{ch,T) - c.{ch,T) (fig. |). 

At the initial fast state of the system relaxation only the phase boundaries attain new equilibrium 
and the defect activation is completed, whereas the state of the phase bulk remains unchanged. The 
latter enables us to adopt the following initial conditions for the Er concentration inside the phase 
Qpd: 

c{z,t) = Co for ;z > 0, t = 0. (6) 

The system of equations (|l])-(P) forms the required relaxation model for Pd-Er-H system. Let 
us discuss its features in detail. 

The proposed model contains four principal parameters: 

^, ^, fi=*^, A = ^. (7) 
Co Co DTtr qoa 

For the sake of simplicity we rewrite system (|T])-(^) in the dimensionless form by introducing the 
normalized variable t] = c/cq and 6 = q/qo, and the dimensionless spatial and temporal coordinates 
e = and r = m/K?'- 

(s) 

where i^tr = Ztr/z^j., and = vz^^/D is the dimensionless velocity of the interface T. The scales of the 
spatial coordinate and time are chosen so that the unit value of ^ correspond to the defect position 
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at initial time and the unit value of r characterize the diffusive penetration of Er atoms over 



scales about z^^.. Equation (H) is subject to the boundary conditions: 



A 



Pd 



lC=o 



Co 



dr] 



A 

Co 



and to the initial condition at r = 0: 



r/(0 = 1 for ^ > . 



(9) 



(10) 



In the chosen frame the "motion" of the trap-defect Ctr(T) and its capacity 6{t) obey the equations: 



tr 



dr 



dO 

-d and — = —hTtOriizt^^ r) . 
dr 



:iv. 



If we ignore the influence of the defect on the interface motion and on the Er atom distribution 
near T (which is well justified at the beginning of the relaxation process for r ^ 1) then equation (|^) 
under conditions (|9|)- ([T0|) admits the automodel solution: 



r/*(e,r) = l-v^-t9oexp (i^l) 
Co 



1 — erf 



where parameter -t^o is the root of the transcendental equation 

Aexp {^D [1-erf (A)] 



Apd 



(12) 



(13) 



and the interface velocity is given by the ratio = ^Q/^JT. Numerical analysis has shown, in 
particular, that for A/cq = 1.5 and Apd/co = 0.5 the value do ~ 0.2 and on time scales t < 1 the 
interface velocity T is also about 1. Therefore, on one side, for the given values of the parameters 
the interface T is to reach the defect in time r being about 10. 

On the other side, if the parameter i7 ^ 1, because the characteristic time rtr of trapping Er atoms 
by the defect is small, and its initial capacity go is such that A < 1 then, by virtue of equation (H), 
the defect will trap all the Er atoms located near it within a sufficiently short time r ~ As the 
result the distribution of Er atoms in the vicinity of the defect will take the form shown in fig. ^ 
and it is characterized by small values of variable ri{ztriT). As follows from the inequality ^ 1 
and equations (|TTp this form does not change until the interface T comes sufficiently close to the 
defect or until the defect is filled up with Er atoms. As it must, the width of the defect neighborhood 
depleted of Er atoms increases in time as ^/T. Therefore, for time scales <^ r ^ 1 and the 
defect neighborhood of thickness about several equation (||) gives the distribution of Er atoms 
approximately of the form: 



^) = erf 



2^7 



and, in addition, the dimensionless rate of Er atom trapping is estimated as 



dine7]5{i-itr)^^v{zt..r)^2 



Qtr 



dr] 



=5tr+0 



vrr 



(14) 



(15) 
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Whence, in particular, setting T]{ztr, r) ^ 1 we have got the estimate of the time during which 
distribution (|14D develops. 

After a lapse of time about t < 1 the regions near the interface T and the defect where the 
Er concentration exhibits substantial spatial variations begin to overlap and the resulting form of 
the Er distribution is shown on fig. It is essential that in the given case the gradient of the 
Er concentration at the interface T changes the sign, causing the interface to move in the opposite 
direction. As time passes further and the defect is filled with Er atoms its effect on distribution of 
Er atoms becomes weaker and the Er concentration in its immediate neighborhood starts to grow 
(fig. ^). Then as follows from the second equation of (pUD and estimate (|l^) after a time r ~ A~^ 



the defect has practically no effect on diffusion of Er atoms and their distribution takes again the 
automodel form (0) (fig. ^) . In this case the gradient of Er concentration at the interface T changes 
the sign again and the interface motion returns to the initial direction. 

The described process is the basis of the proposed model for the nonmonotonic relaxation in 
system Pd-Er-H. So, to substantiate our speculations we present also curves describing the interface 
motion that have been obtained by solving equations (|])-(|ll|) numerically. Fig. ^ illustrates the 
interface motion for various values of the defect capacity in the interval 0.1 < A < 0.8, with the 
other parameters being fixed, A/cq = 1.5, Apd/co = 0.5, and Q = 1000. When the defect capacity 
is large enough so the parameter A < 0.8 we see the nonmonotonic motion of the interface T. 
We have taken the abscissa as the square root of time, a/t, to make the results more clear. It is 
due to the fact that the second stage of the interface motion, i.e., when the interface moves in the 
direction opposite to the initial one, is characterized by the duration exceeding substantially the one 
of the first initial stage. It should be noted similar behavior of the relaxation processes has been 
observed experimentally []TU|, where each of the points of time variations in the phase relative 



volumes changing sign is distant from preceding one by a more and more prolonged time interval. 
We point out once more that obtained results are of qualitative level and to compare quantitatively 
the experimental data and theoretical results a more sophisticated model is required. 
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two phase areas one phase area 




Figure 1: Supposed state diagram for the triple system Pd-Er-H in the region of small Er atom 



concentration. In drawing this diagram we have based on the results of the paper [0 for Pd-Er 
system and on the general assumptions about the hydrogen effect on the state of the Pd-Er alloy. 
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Figure 2: Illustration of Er distribution at the beginning of the relaxation process. 
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Figure 3: Illustration of evolution of Er distribution inside the phase Qpd during the interface motion 

{ta <tb<tc< td). 
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Figure 4: Results of numerical simulation of the interface motion (in dimensionless form): a) the 
interface velocity ^ vs. the square root of time \/r, b) the interface position, {}{T')dT', vs. square 
root of time y/r. (A/cq = 1.5, Apd/co = 0.5, Q = 1000) 
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